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Abstract 

We study the asymptotic limit for large matrix dimension N of the partition 
function of the unitary ensemble (/? = 2) with weight 

/ z 2 t x 2 
w(x) := exp \-— 2 + --- 

We compute the leading order term of the partition function and of the coefficients 
of its Taylor expansion. Our results are valid in the region c±N~2 < z < C2-/V4 . Such 
partition function contains all the information on a new statistics of the eigenvalues 
of matrices in the Gaussian Unitary Ensemble (GUE) that was introduced by Berry 
and Shukla [2|. It can also be interpreted as the moment generating function of the 
singular linear statistics 

N 

E 

3=1 
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1 Introduction 

1.1 Background 

In Random Matrix Theory partition functions of ensembles whose probability measure is 
invariant under conjugation by unitary matrices (unitary or (3 = 2 ensembles) are integrals 
of the form 

— / JJwfo) Yl \%k-Xj\ 2 d N x, (1.1) 

' ^ JN j=l l<j<k<N 

where w(x) > is a weight function and usually J is either an interval, or the whole 
real line or the unit circle. The theory of orthogonal polynomials (see Szego [15], pp. 
23-28) implies that such integrals are determinants of Hankel or Toeplitz matrices, or a 
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linear combination of the two. Computing asymptotic formulae of such determinants is a 
very important task — often a very difficult one — in many branches of mathematics and 
physics. The asymptotics of the integral (11.11) depends crucially on the analytic property 
of w(x). Usually, singular weights are the most challenging. 

The purpose of this article is to compute the following expectation value over the 
Gaussian Unitary Ensemble (GUE): 



E N (z,t) := / [ TTexp 

j* n Vj=i 



+ — ) )P G VE(xu...,x N )d N x, [1:2} 



where 



d gve(xi,...,x n ) : = Y~N\ exp \-\51 x ) ) II \ x k~ x i\ 2 (1-3) 



j=l / l<j<k<N 

is the joint probability density function (j.p.d.f.) of the eigenvalues and 

N 

1 / / 1 ' 

' J N - z 



N-l 



j=l / l<j<k<N 



;i-4) 



is the partition function of the GUE. 
Let 

N 

Ajv(x) := JJ^-Xj-), (1.5) 

3=1 

be a polynomial of degree N whose roots x\, . . . , are all real and define 

Qn(x) " A%(x) -A%)A%(xy (L6) 

This function was studied by Tuck [16] in a numerical investigation of the zeros of the 
Riemann zeta function C( s )- hi Tuck's article A N (x) was replaced by the Hardy function: 

Z(t) := fV'exp (f ) Tr-^r (| + C (| + it) . (1.7) 

The Hardy function is an entire function of order one; it is real for t e R and its zeros 
coincide with non trivial zeros of ({s). The motivation to study (jl.6p was that the Riemann 
hypothesis implies that 

W(t) := Z' 2 {t) - Z(t)Z"(t) > 0. (1.8) 

More generally, if Z(t) is an entire function of order less or equal to one, real on M. and 
whose zeros are all real then W(t) > 0. 
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Suppose that the roots of A^(x) are N random variables with j.p.d.f. p(x\, . . . , xn), 
and denote by P{Qn) the probability density function (p.d.f.) of the random variable 
Qn(x). Berry and Shukla [2] observed that Qn{x) is a sensitive indicator of the degree of 
repulsion between neighbouring zeros of A N (x). More precisely, the rate of decay of P(Qn) 
as Qn — > oo is related to the rigidity of the roots of A N (x). They also studied in detail 
P(Qn) in two cases: when the roots of An(x) are N independent identically distributed 
(i.i.d.) standard normal random variables and when An(x) is the characteristic polynomial 
of a matrix in the GUE. 

Berry and Shukla showed that all the information on P(Qn) is contained in the expec- 
tation value 




in the sense that all the moments of the distribution P(Qn) can be extracted from its 
knowledge. They computed the average (11.91) and the moments of P(Qn) when the roots of 
An(x) are i.i.d. standard normal random variables, but not when they are the eigenvalues 
of matrix in the GUE. 

The integral in the right-hand side of equation (ll.2p is amenable to other interpretations. 
For example, if we set 

s = z 2 /2 = -t, s> 0, (1.10) 

then 

M N (s):= y^JJexp^-s (J^ + ^Jp G ije(xu- ■ ■ ,x N )d N x (1.11) 
is the moment generating function of the p.d.f. of the singular linear statistics 




(1.12) 



Furthermore, Z^E^{^z^€) is the partition function of the unitary ensemble with weight 

Now, denote by TCj(x), j G Z+, the monic polynomials orthogonal with respect to w(x). 
The integral E^(z,t) can be rewritten as a Hankel determinant: 

N-l 

E N (z, t) = Z^ 1 det (p j+k )^k=o = Z n X II h S> ( L14 ) 

3=0 



where 



/oo 
w(x)x^dx, j G Z + 
■oo 



(1.15) 
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and 

/oo 
w(x)TCj(x)TCk(x)dx = hj5jk- (1.16) 
-oo 

If we know the behaviour of the polynomials ttn{%) as N — > oo (Plancherel-Rotach asymp- 
totics), then we can extract information on the asymptotic limit of En{z, t). Our approach 
consists in studying the solution of the Riemann-Hilbert (R-H) problem associated to the 
polynomials 7Tn(x) in the limit as iV —>■ oo. The main tool is the nonlinear steepest descent 
analysis developed by Deift et al [HI E]- The average E N (z,t) can then be computed in 
terms of such a solution using a set of differential identities introduced by Bertola et al [3] . 

After this work was completed, we discovered that independently Chen and Its [5j [6] 
studied the orthogonal polynomial problem associated with the partition function 

i r N 

H N (a,s):=— Yl x l e ~ Xj ~ S/Xj II W~x j \ 2 d N x, (1.17) 



j=l l<j<k<N 



where a > — 1 and s > 0. The weight in this integral is that one of the Laguerre polynomials 
perturbed by the singular factor e~ s / x . In a first paper Chen and Its |5J proved that 
H]y(a,s) can be expressed as the integral of the combination of particular third Painleve 
functions; in a second article [U] they derived asymptotic formulae for the orthogonal 
polynomials, the corresponding recurrence coefficients and the h^s. When the parameters 
t and a in equation ( 11. 131) and x a e~ x ~ s / x respectively are both set equal to zero, then the 
systems of orthogonal polynomials associated with the two weights can be mapped into 
each other by the change of variables x t— > x 2 /2. The Plancherel-Rotach asymptotics of the 
orthogonal polynomials can be studied using the nonlinear steepest descent in both cases. 
However, even when t — a — 0, the partition functions (II. 2p and (11 . 17)) are not equivalent 
and cannot be mapped into each other by a simple change of variables. 

1.2 Statement of results 

The average ( 11.2)) is an entire function of t, thus its Taylor series has an infinite radius of 
convergence and we can write 

oo 

E N ( Z ,t) = J2 E Nm(z)t m . (1.18) 

The main goal of this paper is to compute the leading order term in the asymptotic ex- 
pansion of Etf(z,t) and EN m (z) as iV — ► oo. Our main result is the following: 

Theorem 1.1. Let C\N~^ < z < C2N^ , where C\ and Ci are two constants independent of 
z and N. The expectation value E N (z,t) is given by 

, s ( z 2 9 / 2 4 \ t 2 Nl 

E nM = Bpj exp ^— - -y (N^ - l) + (i u) 

x (l + o(l)), iV^oo, 
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where Bn is the ensemble average 

N 



B n-= J rN (JI e ^JP G VE(x 1 ,...,x N )d N x, (1.20) 

which is independent oft. 

Berry and Shukla [2] showed that the m-th moment of P(Qn) is given by the integral 

2m /»oo 

M Nm := 1 x ~ m Y[nl z 2m ~ x E N2rn {z)dz. (1.21) 

n=m ^ 

From Theorem 11.11 it is straightforward to compute the coefficient Ex2m{z) in the series 
expansion (11.181) . 

Corollary 1.1. Let c\N~^ < z < c 2 Ni , where c\ and c 2 are two constants independent of 
z and N. The leading order term of the coefficient of t 2m in equation (11.181) is 

Z -7-^w{ Niz ^ ~ l ) ) ] ^ N^oo. (1.22) 

4 2 3 v V 2 3 mlz 3 

Unfortunately, the asymptotic limit in equations f l 1 . 1 9 f) and ( 11.221) cannot be assumed 
to be uniform in z: there may be non negligible contributions from the region z < CiN^z 
that would affect the integral (1 1 . 2 1 D . Such contributions can be investigated by studying 
the double scaling limit of the matrix ensemble with weight ( 11.131) . This will be the subject 
of a forthcoming publication. 

The structure of this article is the following: in £0 we introduce the R-H problem for 
the orthogonal polynomials with weight (11.131) (after appropriate rescaling of x, z and t) 
and the differential identities used to compute the leading order asymptotics of E N (z,t); 
in ^H] we find the equilibrium measure on which the R-H analysis is based; in §3] we apply 
the nonlinear steepest descent to the R-H problem; §5] and £0 are devoted to complete the 
proof of Theorem 11.11 combining the asymptotics of the solution of the R-H problem and 
the differential identities discussed in 32l 
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2 Preliminaries 

Let us introduce the scaling 



(|) 2 and Vs = ^L, j = l,...N. (2.1) 



t f "\ i i 

Vl = — F=i 1>2 



The weight ( I1.13P becomes 



^):=exp(-iV(0 + |)+|). (2.2) 



We also define the partition function 

1 „ AT / , 2 \ \ AT 

which is proportional to the average E^(z J t) J namely 

E N (z,t) = Z N 1 N-G N (v 1 ,v 2 ). (2.4) 

For convenience, where there is no risk of confusion with the quantities introduced in §l.ll 
we denote the polynomials orthogonal with respect to wn{v) by itj(y). Similarly, we write 

N-l 

Gn^, v 2 ) = det (fJ, j+k )J~^ = Y[ h v ( 2 - 5 ) 

where 

/CO 
w N (y)y 3 dx, j G Z+ (2.6) 
-co 

and 

/•CO 

WN(y)^j(y)^k(y)dy = hj5 jk . (2.7) 



Let us define the matrix valued function 



y(v] ■= I ™ 2 ™ J -°° s ~y I (9 ? 

Vi/J • I / \ WjV-l TOO 7T J v_ 1 (s)u. J y(s) ,^ I 5 

^-l^tf-l 12/ j ^— J.^ — dS y 

where «A7_i = —li^ih^-x. Fokas et al pUl EH] showed that ^(y) solves the following R-H 
problem: 

1. Y(y) is analytic in C/R, 

2. = (J «"f), tf6 R, (2 9) 

3. Y(y)=(I + 0(y- 1 ))^ y _ oo 
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where Y + (y) and Y-(y) denotes the limiting values of Y(y) as it approaches the left and 
right-hand side of the real axis. It turns out that the partition function Gn(vi,v 2 ) can be 
expressed in terms of Y(y). 

Lemma 2.1 (Bertola, Eynard and Hanard |3J). The following differential identities hold: 

d log G N 



dvi Airi 
d\ogG N N 



-I - Tr (Y-\y)Y'(y)a 3 ) dy, (2.10a) 

J y =o y 

<f \Tr(Y-\y)Y'(y)a 3 )dy, (2.10b) 
Jy=0 y 



dv2 Swi 

where the contour of integration is a small loop around y = oriented counter-clockwise. 

Proof. Taking the logaritmic derivatives of both sides of equation (12.51) and using the 
orthogonality conditions (12 .7p gives 



3=0 " JJ 



These integrals can be rewritten as 



dlogG N [°° K N {y,y) 

-^r =[ ~ N) L^^ d% k=1 ' 2 ' (2 - 12) 



where K^(x,y) is the kernel 



K N {x, y) := Vw N (x)w N (y)J^ l^llll (2.13) 

3=0 3 



N-l 



and 

G N {v u v 2 ) = \ N\ J| h d J det(K N ( yj ,y h ))"-' . (2.14) 
In order to evaluate the integral fl 2 . 1 2 [) we use the relation 



y/w N (x)w N (y) ( x fQ 1 



Tr( Y~ 1 (y)Y + (x) ([ 
2m(x -y) \ + yy ' +K J \0 



(2.15) 



which follows from the Christoffel-Darboux formula, the definition of Y(x) and the fact 
that dety(y) = 1. By using l'Hospital's rule we obtain 

K N (y,y) = ^Tr(Y-\y)Yi(y) (° ^ ). (2.16) 
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Combining the jump condition of the R-H problem (12.91) and equation (12.16j) gives 

K N (y,y) = ^(Tr (Y- 1 (y)Y^y)a 3 )) - Tr {Y-\y)Y[(y)a,)^ y eR. (2.17) 

The asymptotic behaviour of Y(y) as y — > oo implies that the functions 

Tr (Y ± l (y)Y^y)a 3 ) (2.18) 

are analytic in the upper/lower half planes with a simple pole at infinity. Finally, the 
identities (12.101) follow from the residue theorem. □ 

This lemma gives an explicit link between the solution of the R-H problem and the 
average (11.21) . The main challenge that we are facing is to compute an asymptotic formula 
for Y(y) using the nonlinear steepest descent method. Then, we can obtain a formula for 
E N (z,t) using equations (12.101) . 



3 The equilibrium measure 

The asymptotics analysis of Y(y) relies on computing the g-function and the support of 
the equilibrium measure for the potential 

where we have neglected the term v±/y because it is asymptotically small. The equilibrium 
measure is the positive normalized Borel measure fi(y) that minimizes the energy function 



OO POO 



H^) : =- / log \x - y\dfx(x)dfi(y) + V (y)du,(y). (3.2) 

J —oo J —oo J —oo 

It satisfies the conditions 

/oo 
log \y - s\dfi(s) - V (y) = I, x E Supp (/i) , (3.3a) 
■oo 

/oo 
\og\y - s\dfi(s) - V (y) < I, i6l/ Supp (//) . (3.3b) 
■oo 

for some constant I. Moreover, if fj,(y) satisfies (13.31) . then it must be equal to the equilib- 
rium measure. 

The g-function is defined by 

POO 

g(y):= / log(y - s)dfi(s). (3.4) 
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It is analytic in C \ R and has jump discontinuities on the real axis. The conditions ( 13. 3ft 
expressed in terms of g(y) become 

9+(y) + 9-(y) - V (y) =1, ye Supp (ji) , (3.5a) 

Re(g+(y) + g^{y)) - V (y) <l, y E R/ Supp (/x) . (3.5b) 

Furthermore, we have 

g(y) = logy + Oiy^ 1 ), y -> oa. (3.6) 
We shall derive an expression for g(y) using the ansatz 

5 U/J - -3- ^ • ^ 
The constants Ai, A2 and A3 are determined by the asymptotic behaviour of g'(y), namely 

g'(y) = - + 0(y- 2 ), y^oo, (3.8a) 

y 

g'(y) = 0(l), y^O. (3.8b) 

These constraints give 

X\ + \ (^2 + Ajj) = 2, (3.9a) 

Ar 2 + \ (A2 2 + A3 2 ) = 0, (3.9b) 

A?A 2 A 3 = -ua. (3.9c) 

Let Aj = A|. Then, we see that A 2 , A 3 are the solutions of the quadratic equation 

y 2 - 2(2 -A 1 )y + ^ = 0, (3.10) 



while v4 2 1 and A 3 1 are solutions of 



2 



2 A* 



^ + -^ + -^ = 0. (3.11) 

From equation (13.101) it follows that 1 and A^ 1 also satisfy 

^-2(2-^ + 1 = 0. (3.12) 

By comparing this equation with ( 13.111) . we see that At is a solution of the equation 

A\-2A\-vl = Q. (3.13) 
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Lemma 3.1. There exists a solution to equation (13. 9ft such that Ai G iR and A2, A3 G 
Proof. The solutions of (13.111) are 




(3.14) 



Therefore, if a solution of (I3.13P such that — < < exists, then the two solutions 

to (I3.14p are real and positive and the lemma is proven. Now, if A\ = 0, then the left-hand 

3 

side of (I3.13P is —v\ < 0; if A± — — it is 2t>| > 0. Hence, there is a solution to ( I3.13P 
between and 0. □ 



Let us now choose Ai, A2 and A3 such that Im(Ai) > and < A2 < A3. Furthermore, 
define 

)H: ^_^t«pF3, + i ( , 15) 

where the integration path is chosen such that it does not intersect the interval (—00, A3) 
and 



, = - 2 . im 1 m - logy - r c - MK/g - WEEM ds 1 . ( 3. 16 ) 



The function g(y) is analytic in C \ (— 00, A 3 ) and g(y) ~ logy as y — ► 00. Now, let 
£ := Si U £2, where 

Ei := [-A 3 , -A 2 ] and E 2 := [A 2 , A 3 ] . (3.17) 

We will dedicate the rest of this section to show that g(y) satisfies the conditions (13. 5p 
with 

Supp (//) = E = Ei U E 2 . (3.18) 
Lemma 3.2. The function g(y) defined in equation ( 13. 15ft satisfies the jump discontinuities 

9+(v)+9-(v)-V (y) = l, 2/GE, (3.19a) 
9+(y) ~ 9-(y) = 2n h ye(-oo,-A 3 ), (3.19b) 

9+(y) - 9-{y) = ni, y€(-A 2 ,A a ). (3.19c) 

Proof. Since Vo(y)/2 has no jump discontinuities in C, we only concentrate on the integral 
in the right-hand side of (13.151) . Let us write 



(y 2 -A?)vV^Ai)(^-AL 
2y3 



u(y) := ^ CUXML C*L, ( 3 . 2 0a) 



g(y) := I v{s)ds. (3.20b) 

JX 3 
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If y G S 2 , then 

9+(y)+g-(y)= I +f-00) rfs. (3.21; 



a 3 

Since z/(y) changes sign across S, we have 

v+(y) = -v4y), yGS. (3.22) 

Equations (13.211) and f)3.22p imply that (I3.19ap is satisfied for y G S 2 . 

Suppose now that y G Si. Let T± be a contour that consists of 3 parts: the first one 
goes from A3 to A2 on the positive/negative side of the real axis; the second part T± is 
a semicircle from A2 to — A 2 in the upper/lower half plane; the last part T± goes from — A2 
to y on the positive /negative side of the real axis. Then, the jump on Si is given by 



9+(y) + 9-(y) = / v(s)ds+ / v(s)ds= / (f+(s) + v-{s))ds 
Jr+ Jr. J-\ 2 

+ / (v+(s) + v_(s))ds + / v(s)ds+ / u(s)ds. 
J\ 3 Jr 2 + Jr 2 _ 

Since v+{s) = —V-(s) on S, the first two terms in the right-hand side of equation (I3.23f) 
are zero. Furthermore, under the map s t— > — s, the contour r+ becomes — and v(—s) = 
—v(s). Hence, the sum of the last two terms are zero too. This proves (I3.19ap on Si. 
Let y G (—00, — A3). We simply have 

g+(y)-g-(y) = 2mHeag(y). (3.24) 



y=co 



Then, equation ( 13.19bl) follows from 

*»M = - X ' + M + ® ,2 = -l, (3-25) 

y=co Z 

where we have used (13.9aj) . 

Let y G (— A 2 , A 2 ) and consider the contours T± := T± U T± and F±. Now r± joins A 2 
to y and and r± connects — A 2 to —A3 respectively. As before ± denotes the upper/lower 
half plane respectively. We have 

g+(y) ~ g-(y) = / v(s)ds- \ v{s)da = \ v(s)ds, (3.26) 

Jf+ Jf- Jr\u(-rl) 

where we have used 

_ . . A7 2 + (Xo 2 + \^ 2 ) 12 . . 
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which follows from equation (13.9bl) . The function u(s) is odd and T± is mapped into — 
by s t— > —s, therefore 



u(s)ds= / u(s)ds. (3.28) 

f3_u(-fi) A^u(-r^) 

Finally, Cauchy's theorem gives 

/ v{s)ds = iri Res g(y) = —in. (3.29) 

^u(-rL) y=°° 

This completes the proof of equation ( I3.19cl) . □ 

We are now ready to prove that the definition (13. 15ft gives the g-function. 

Proposition 3.1. Suppose that Ai, A2 and A3 satisfy the conditions of Lemma \'S.1\ that 
Im(Ai) > and < A2 < A3. Then, the function g(y) defined in ( 13. 15ft satisfies 

g+(y) + g-(y)-V (y) = l, y e E, (3.30a) 

Re (g + (y) + g.(y)) - V Q (y) < I, y G R \ S. (3.30b) 

Proof. Equation ( I3.30al) was proven in Lemma I3T21 We are left to prove inequality ( 13.30bl ). 
From the jump discontinuities (I3.19p . we see that outside S the real parts of g+(y) and 
g~(y) are equal. In particular, the real part of g(y) is continuous outside S. 
Equation (I3.19al) and 

<Kv) = ^-s(y) + 5 (3-31) 



give the equality 
Now note that 



Be(g+(y)+g-(y)) = Q, V e S. (3.32) 

(3.33) 



u(y)>0 if y E (-A 2 ,0)U(A 3 ,oo), 
u(y)<0 if y G (-00, -A 3 ) U (0, A 2 



Thus, Re^(y) is an increasing function in (— A 2 , 0) U (A3, 00) and decreases in 
(—00, — A3) U (0, A 2 ). This property and (I3.32f) imply 

Re(g+{y) + g^y)) > 0, xGl\E. (3.34) 

The proposition now follows from equations (13.311) and (I3.34H . □ 



4 The Riemann-Hilbert analysis 

The purpose of this section is to study the asymptotic limit of the solution Y{y) of the R-H 
problem ( 12. 9ft using the nonlinear steepest descent analysis developed by Deift et al [H [9] . 
One of the main ingredients is the g-function computed in §HJJ 
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4.1 Deformation of the Riemann-Hilbert problem 

Let us introduce 

F(V) := (i sq + (s)(s-y) + * / A2 ' (4,1) 

where q(y) and £ are denned by 

?(!/) == )/(f 2 - A i)(f 2 - A l)' ( 42a ) 
K := 2 / --, (4.2b) 

_ / s ^y 2™_ 

J-\2 q{s) 

The right-hand side of equation (I4.2cl) follows from the residue theorem. 

The function F(y) is bounded at the points ±A2, ±A 3 and satisfies the following scalar 
R-H problem: 

1. F(y) is analytic in C \ [—A3, A3], 

2. F+ (y) + F4y) = ^, y E S, 

y (4.3) 

3. F + (y) - F4y) = £ Vl , ye(-X 2 ,X 2 ), 

4. F{y) = 0{y- 1 ), y -> 00. 

Let us define 

T(y) := e (-^)- 3 y (y)e -(iV S ( 2/ )-F( 2/ )).3 e ^ ; (4 4) 

where = (0 -1 )■ The matrix function T(y) is the solution to the R-H problem 

1. T(y) is analytic in C \ R 

2. T + {y)=T„(y)J T {y), y e R, (4.5) 

3. T( 2/ ) = / + 0( 2 /- 1 ), l/^oo, 

where 



-iv( s+ (y)- 5 _(y))+F+(y)--F-(!/) p -N{g+{y)+9-{y))-F+(y)-F-{y)+^ , 
q e ^^+(w-s-(wJ--p+(w+- F -(w 



and is defined in (13. 20b j) . 



We now perform a standard technique in the steepest decent analysis (see pfl,[8j[9]): we 
open lenses around the intervals Sy, j = 1,2. The interiors L±j and contours ELy of the 
lenses are defined as in figure HJ 
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Figure V. The opening of the lenses around the intervals £j, j = 1,2. L±j and H-y are 
the interiors and contours of the lenses in the upper /lower half plane respectively. We also 
denote Ej = H + j U H_j and H = Hi U H 2 , j = 1, 2. 



Let us introduce the matrix function 
(T(y), 



S(y) ■= 




_ e 2Ng(y)+2F(y)-^ } 



e 2Ng(y)+2F{y)-^- - { 



for j = 1,2. It satisfies the R-H problem 



1. is analytic in C \ 

2. 5 + (y) = 5_(y)J 5 (y), 

3. S(y) = I + 0(y' 1 ), 



where Js{y) is the jump matrix 



My) < 



yeC\ (L +j U L_,-) , 



2/ -> oo. 



e 27V9(j/)+2Ffe)-^ x 

i N 

-1 

'j e -2iVs(2/)-2F(2/)+a- > 

1 

,i\M+£t)i e -A r (s+fe)+9-(j/))-^+(?/)-^-(?/)+^ 

Q e -NTTt-(vi 



y e H, 
yes, 

x G R \ (-A 3 , As) 
y G (-A 2 , A 2 ) 



(4.7) 



(4.8) 



(4.9) 



The conditions (13.301) imply that away from some small discs D±\ j of radius 5 centered at 
±Aj, j = 2, 3, the off-diagonal entries of the jump matrix Js(y) are exponentially small in 
N, except on the intervals Si and £ 2 . This suggests the following approximation to the 
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R-H problem for S(y): 

1. S°°(y) is analytic in C \ [—A3, A3] 

2. S?(y) = S?(y)J">{v), 

3. S°°(y) = I + 0(y- 1 ), 



where 







J°°(y) ■= { 



1 

-1 0. 



y e [-A 3 , A 3 ] 
y^oo, 

y E (-A 2 , A 2 ), 



(4.10) 



(4.11; 



The approximation S°° (y) is known as outer parametrix. 
4.2 The outer parametrix 

Here and in the rest of §U ±\j will always refer to the edge points of the support S of the 
equilibrium measure. Hence, j — 2, 3 only. 

The solution to the R-H problem (14.101) exists and is uniformly bounded in N outside 
of small discs D±\. around the points ±A 3 -. Such a solution can be constructed in terms of 
elliptic theta functions as in Deift et al [9]. Here we follow their treatment. 

Let £ be the elliptic curve 

q 2 = {y 2 -\l){y 2 -\l) (4.12) 

and choose a canonical basis of cycles as in figure EJ Then, the holomorphic 1-form oo{y) 

a 

b 




Figure 2: The a and b cycle of the elliptic curve (14.121) . 
dual to this set of cycles is given by 

/ \ dy 

K q{y) 

where Kq was introduced in equation (14.2bjl . The Abel map is defined by 

" y ds 



u(y) 



a 3 



K q(s) 



(4.13) 



(4.14) 
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where the contour of integration is chosen such that it does not intersect the interval 
(— oo, A3). Let II be the 6-period of the 1-form ui{y): 

n : = 2 dy / - . (4.15) 

The elliptic theta function for the curve (14.121) with this choice of cycles is given by 

9(s) := ^ e wIIm2+2mm , (4.16) 

Consider the function 



7 ■= (i^k±M ) 3 (4 it) 

where the arguments of the individual factors in the fourth root are chosen to be between 
— 7r and 7r. Then the solution to (I4.10p is given by 



7+7" 1 9 




2/rl 


+d) 


2 

7 _ 7 -i e 


<9(tt(yH 
(«fe)-f- 


-d) 
"1? 

2iri 


~") 


2i 


e(«d/)- 
0(u(oo) 


-d) 
1 + 


d) 




:oo)-f 




2wi 



Coo/ A — o-l 2 -2i 8(-u(y)+d) I fA 1S fl \ 

2 0(u(,y)+d)* 



H := diag ( ? rA^T^v fl(tt(oo) + d) ( 

1 x !f + d) , e(-u(oo)-f-^ + d) 1 



where d is the constant 



in . , 1 n f X2 ds f° ds , A , 

2 2 2 2 J Xa K q + {s) J x , 2 K q + {s) 

Using the definition of 6-period (14.151) . of K (I4.2bl) and the fact that q(—s) = q(s) give 
d = -\. 



4.3 Local parametrices near ±A2 and 1LA3 

Near the edge points ±Aj the approximation of S(y) by S°°(y) fails. Therefore, we must 
solve the R-H problem (14. 8 j) in small neighborhoods of these points and match the solutions 
to the outer parametrix (14.101) up to an error term of order 0(A r_1 ). More precisely, let 
5 > and D±\. be a disc of radius 5 centered at ±A 3 -. We would like to construct local 
parametrices S( ±Xj ^(y) in D±\ j such that 

1. S^'^) is analytic in D ±X] \ (D± Xj fl (R U S)) , 

2. S^\y) = S {± ^\z)J s (y), y G D ±A . n (R U S) , (4.20) 

3. S^y) = (J + O^- 1 )) S°°(y), y G <9D ±A .. 
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These local parametrices are given by 

S^\y) := E^\y)¥^\C±j)e^ Nm+F(v) -^. (4.21) 

The matrix functions \I/( ±a j) are constructed using Airy functions (see, e.g., [7] pp. 213 
216). The explicit expression of \I/( ±A ^ is quite lengthy. Besides, it does not enter in our 
calculations. Therefore, we refer the interested reader to the original literature (see [H El 
EJ[9]). The C±Aj's are conformal maps inside the neighborhoods D±\ j given by the analytic 
continuation of the functions (± x to the whole D±\. : 

Ctx 3 ■■= (l^y (g(y) ~ ~g(X ±J ))K lm(y) > 0. (4.22) 

The matrix En S \y) is invertible and is the analytic continuation to the whole D±\ of 
the following quantity: 

E^\v)--=^e^S-{y)e^^ (±( ±Xj f. (4.23) 

Remark 1. Since E^t (y) is analytic inside D±\ ., we see that near the points ±Aj, the 
function S°°(y) behaves like 



S^(y)^S^\y)(yTX j )-^ 



1 ±(-l)^+ 1N 

T(-iy +i i ; (4.24) 

Xe -^ e (-%)+|)^ y^±X j , 



where SQ 3 \y) is holomorphic and invertible at ±Aj. 

4.4 The final transformation of the Riemann-Hilbert problem 

We now show that the parametrices we constructed in §4.21 and §4.31 are indeed good 
approximations to the solution S(y) of the R-H problem (14.81) . 
Let us define 

m:=[ s ^ s ^W l - veD±x - (4-25) 

m \S(y) (S"(y))' 1 , yeD ±Xj . 1 ; 

Then the function R(y) has jump discontinuities on the contour Tr shown in figure [SJ In 
particular, R(y) satisfies the R-H problem 

1. R(y) is analytic in C \ Tr, 

2. R + (y) = R_(y)J R (y), yeT R} (4.26) 

3. R(y) = I + 0(y~ 1 ), y - oo. 
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Figure 3: The contour Yr. 



From the definition of R(y) it follows that the jump matrix Jr{d) has the following order 
of magnitude as iV — > oo: 



Mv) 



I + OiN- 1 ), yedD ±Xj , 

I + O (e~ Nri ) , for some fixed rj > for y 6 T R \ D±\ j 



(4.27) 



Then, using well established techniques (see,e.<7 [8] §7) we obtain 

R(y) = I + ' ' 



N(\y\ + 1) 



(4.28) 



uniformly in C. Therefore, the solution S(y ) of the R-H problem H4.8[) can be approximated 
by S°°(y) and S^(y): 



S(y) 



(i + oiN-^s^iv), y^D ±X]1 

(I + OiN-^S^y), yeC\D ±Xr 



(4.29) 



Combining these expressions with equations (14.71) and (14. 4p we obtain an asymptotic for- 
mula for the solution of the original R-H problem H 2 . 9 f) . which can be inserted in the 
differential identities (12.101) . 



5 Asymptotics of the differential identities (12.101) 



In this section we will compute an asymptotic formula for the partition function G^(vi, t^) 
defined in ( 12. 3ft . The analysis is similar to those ones carried out in [12, 



5.1 The finite v<i regime 

Consider the trace in the integrals (12.101) : 

iy):=Tr(Y- 1 (y)Y\y)a 3 ). (5.1) 



Since a(y) is analytic in a neighbourhood of the origin, in order to compute the inte- 
grals (I2.10p we only need the first two terms in its Taylor expansion at y — 0. We have 

a(y) = Tr (S- 1 (y)S'(y)a 3 ) + 2Ng\y) - F\y). (5.2) 
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By using S(y) = R(y)S°°(y) and the expression (14.181) for S°°(y) we obtain 

a(y) = Tr (&' 1 (y)&'(y)a 3 ) + 2Ng'(y) - F'(y) + O f , (5.3) 

where 

Q(y)=H- 1 S°°(y). (5.4) 

It is convenient to split the computation of the right-hand side of this equation in two 
parts. We first determine 

a (y):=Tr(e- 1 (y)e'(y)a 3 ), (5.5) 

then we evaluate 

2Ng\y)-F'{y). (5.6) 

Lemma 5.1. We have 

M . . 

«o(y) = ^y, (5-7) 

where q(y) is defined in (14.121) and M is a constant. 

Proof. We first show that cto(y) is analytic in C \ S and that 

a ,+(y) = -<xo,-(y), y e s. (5.8) 

By differentiating the jump conditions in (14. lip , we see that Q'(y) has the same jumps 
as Q(y). Therefore, the quantity 6 _1 (y)Q' \y) has the discontinuity 

(e- 1 (y)e'(y)) + =(J 00 (y)y 1 (e-\ y )e'(y))_J 00 (y), ye(-\ 3 ,\s) (5.9) 

and is analytic elsewhere. The jump discontinuities of J°°(y) in (14.111) imply equation (15. 8ft 
and that 

a 0:+ (x) = a Q -(y), y G (-A 2 ,A 2 ) (5.10) 

Therefore, a (y) must be equal to q(y) multiplied by a rational function r(y). Since a (x) 
is of order O (y~ 2 ) as y — > oo and is analytic away from the points ±A 2 and ±A 3 , r{y) can 
only have poles at the edge points ±Aj, j = 2,3. From the behaviour of S°°(y) near ±Xj 
(see equation (14. 24ft ). we see that ao(y) can at worst have a singularity of order 1 at ±Aj. 
This implies that r(y) is analytic at the points ±Aj. Hence, it must be a constant M. This 
proves the lemma. □ 

The constant M in (15.71) can be found from the asymptotic expansion of Q(y) as y —>■ oo. 

Lemma 5.2. The function a^{y) in (15. 7p is given by 
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Proof. From equations ( 14. 10ft and (15. 4p the asymptotic expansion Q(y) at y = oo has the 
form 

Q(y) = H- 1 + ^ + 0(y- 2 ), y^oo, (5.12) 

y 

where H is the constant in ( 14.181) . Hence, 

a (y) = -Tr(^± ( r^+O(y- 3 ), y - oo. (5.13) 
The diagonal entries of 0i are 

(0l)22 " "^T [l^MTfTW^) ~ *(«(«>) -i) J +7li " 22 ' (5 - 14b) 

where 71 is the coefficient of y~ l in the expansion of (7 + 7 -1 ) /2. Finally, equation ( 15. lip 

follows by substituting (15.141) into (15.131) . □ 

We can now compute the logarithmic derivatives (12.101) . 

Proposition 5.1. Let he of order 0(1) as N — > 00 and define 

C := ^ r * (5 , 5) 



Then, 

d log Gat 7rzAf 9 



(9f 1 ^2^0^ cfoi 



tog(fi(tt(oo) + ? -^)e(tt(oo)+f + ^) 



_ 2lf "2£)^ + 0(iV ~ 1) ' JV, ~ >0 °' (5 ' 16a) 



lslogG ~ w (i-S(&-i) ,+ a) +0 ^>- (516b) 



N dv 2 



Proof. We are only left to determine the term (I5.6p . First note that the derivative of F(y) 
satisfies 

1. F'(y) is analytic in C \ E, 

2. F;(y) + F:(y) = -^, y E E, ( 5 .17) 

3. = 0(r 2 ), 2/^oo. 
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Moreover, since F(y) is bounded at ±Xj, j = 2, 3, F'(y) cannot have a singularity of order 
higher than 1 at these points. A function with these properties must be of the form 

F'(y) = -^ " lA2A3 ~ , (5.18) 

for some constant C . This means that 

m ■= ; iAaA3 " ^ dy (5.19) 

2|/ 2 vV-A!)(2/ 2 -A!) 

is a meromorphic 1-form with a singularity at the origin such that Q(y) ~ v\dy/2y 2 as 
y — > and is holomorphic elsewhere. To determine the constant C, note that the jump 
conditions (14.31) implies that the a-period of Q(y) must vanish. This gives C = C, where 
C is defined in equation (I5.15p . 

Therefore near y = the function ot{y) behaves as 




+ 0(„ 2 ) + OiN- 1 ). 

Finally, inserting the right-hand side into equations (12.101) gives formulae (15.161) . □ 
5.2 The small V2 regime 

When v% becomes small, the local parametrices constructed in §4.31 must be defined in 
neighborhoods D±\ 2 whose radii 5 is smaller than A2. This affects the magnitude of the 
error term in the asymptotic formulae (I5.16p . We now study its effect. We shall see that 
by constructing local parametrices in shrinking neighborhoods of the points ±A2, we can 
extend the validity of formulae (15.161) for v 2 > N~ 3+e . 
Since the error term is of order 

O I 15 {y)l 6 3 I (5.21) 



IC 



±2 



2 



on the boundary of D±\ 2 , we need to know the order of magnitude of S°°(y) and of the 
conformal map C±A 2 - 

Proposition 5.2. The following asymptotic formulae as v 2 — > hold. 
1. The orders of magnitude of points Ai , A2 and A3 are 

\ 1 = (4) , A 2 = O (4) and A 3 = O (1) . (5.22a) 
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2. Let D±x 2 be small discs centered at ±A 2 with radii 5 < A 2 . Then, on the boundary of 
D±\ 2 the conformal map (±\ 2 defined in ( 14. 22ft is of order 

C±a 2 = O (iVit;~*<y) . (5.22b) 

5. There exists a constant k = 0(1) as u 2 — > such that the outer parametrix S°°(y) 
defined in equations ( 14. 18ft 0/ order 



= O 



X2 
S 



4 fcSL 



e ^ . (5.22c) 



^. T/ie order of magnitude of the exponential in equation (15.211) zs 

exp (2 \F(y)\ + = O (e fc ^) . (5.22d) 

Proof. Equation ( I5.22al) is an immediate consequence of formulae (13. 13f) , ( 13.111) and ( 13.141) . 
As V2 — > we obtain 

Ai = (-2)-*4+0(O, (5.23a) 
A 2 = f~V Ai + Ofa), (5.23b) 



2, 

A 3 = 2 + O (uf) • (5.23c) 
Next, consider the conformal maps C±a 2 - Inside D±\ 2 they behave as follows: 

C ± a 2 = iV§ (cp ± (y - A 2 ) + O ((y - A 2 ) 2 )) , (5.24) 

where 



(AI- A?)t(A 2 (AI-A§) 



<P± '■= ± ^ — ■ (5-25) 

Hence, (±a 2 is of order 



2f A| 



.2 — . 



C±a 2 = iVVi (5.26) 



on the boundary of -D±a 2 - 

Proving ( 15.22c[) requires more work. Firstly, consider 7±7 _1 in a small neighbourhood 
of A 2 . If \y — A 2 | = S is smaller than A 2 , we have 



2A 2 J \A 2 — A3 y \ 8 J V A 2 + A3 / / \ \|A 



(5.27) 
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The case when y = — A 2 can be treated analogously. 

Let us now consider the theta functions that enter in the definition (14. 18ft . First note 
that in the limit as A2 — > 0, the holomorphic 1-form u(y) in (14. 131) becomes a meromorphic 
1-form with a simple pole at y — with residue l/(2ni) (see [JJ and [Hj). More explicitly 

\ 3 dy 



u{y)^- ^==, A 2 ^0. (5.28 

Furthermore, the constant K has the limit 

2; 
A 3 



Ko = -7- + O (A 2 ) , A 2 -> 0. (5.29) 



By writing the period II as 




4 K q + (s) 



(1 + (A 2 )) + /" 2 - - ^== (1 + (A 2 )) 
A? f^AWs 2 - A 2 . 



2 

A3 ./A2 



Using this expression we obtain 



2 ivr \2A 2 / V \A 2 



(5.30) 



'a, ^oSa/s 2 - A| A| ^Asa/s 2 - A2 

we can compute its limit: 

n= (— log |A 2 | -— log|16A 2 |Vl + 0(A 2 )), A 2 ^0. (5.31) 

y 7TI m J 

The Abel map u(y) becomes 

rM ry Yl r y 

u{y) = / u(s) + / u(s) = -+ u(s). (5.32) 



(5.33) 



We can now substitute equations (I5.3ip and ( I5.33f) into the theta functions in (14. 18f) to 
obtain their orders of magnitude as w 2 — > 0. Let A be a constant vector that is independent 
of y and t> 2 . We have 

9(s) = 6 (u(y) + A) = J2 e mUm2+2 ™ {u{v)+A)m . (5.34) 
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The arguments of the exponentials become 

m 2 n + 2 (u(x) +A)m = m(m + 1)11 + 2 [ — ( ) film 



7ri \ 2A 



\2 



3 - 

5 \ 2 



(5.35) 



+ 0( (-) )+0(A 2 ), A 2 ^0. 



The asymptotic behaviour of the period IT in (I5.3ip gives 

9 (u(y) +A) = 1 + e-(^y- 2mA (l + O (J-^j + 0(A 2 )) , A 2 -> 0. (5.36) 
By substituting 

into (15.36p . we see that the matrix elements in ( 14. 18 aft are bounded by infinity and zero. 
Note that, although the term vi^/(2m) depends on t> 2 through £, since our goal is to study 
the coefficients of v±, we can always let v\ be arbitrarily small so that the term vi^/(2iri) 
will only introduce a negligible error into (I5.36p . In particular, we have 

N «i£ 1" 



H<y)-Y-m-l)=0{e^), (5.38) 

where k\ is of order 0(1) in t> 2 . 

From equation ( I5.28P and the fact that A3 = 0(1) it follows that the term u(oo) in the 
constant H in (14.18bl) remains finite as t> 2 — > 0. In fact, we have 

u(oo)=/ [l + O (%))=- + Oft), A 2 -0. (5.39) 

Proceeding as in the derivation of (15.361) . we see that also H remains finite and non-zero 
as f 2 — > 0. 

Arguments similar to those ones that led to equation (I5.33P give 

\F(y)\ = o(h2^ and J = °(hj^, (5-40) 

where y = ±A 2 + 5, and fc 2 and k^ are of order 0(1) in t> 2 . Let us write 

k := max {k\, /c 2 , k 3 } . (5-41) 

Equation ( I5.22cl) follows by combining formulae (15371) . (ESHD and (15391) . and equation (l5~22dl) 
is a simple consequence of (I5.40j) . □ 

We are now ready to prove 
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Corollary 5.1. Let < e < 3. Then, formulae (15.161) hold for v 2 > N 3+e with an error 
terms of order 

£ = ^N~h k ^ . (5.42) 
Proof. Combining equation (I5.2ip with the asymptotic formulae (15.221) gives 

*-°(4£)- <- 3 » 

Then, equation (15.421) follows by setting t> 2 = 0(N~ 3+t ) and 

5 = [N~ l+2 ^ = O (A 2 iV~i) . (5.44) 

□ 



6 Asymptotics of the ensemble average E^(z, t) 



We are now in a position to give an asymptotic formula for the ensemble average (JH 
and complete the proof of Theorem 11.11 

Let us translate formulae (15.161) back into the original variables z = Ny/v 2 and 
t = y/Nvi. We obtain 

aiog G „ z (^_^((^_^Y + ^ ]]+0 ( Z N-^), (6.1a) 

'i 



dz I 2 16 \ \X 2 2 XIJ \j 

d\ogG N iriXf d 



log (e{u(oc) + <;-^ m )9(u(oo)+,+ 



dt ~ v 2 K £ dt 1W& V v v ; T s war; w V" v ^ ^ s ^ wf 

( +ofiV-3-SV (6.1b) 



2iVAf 2A 2 A 3 ViV 

where the Aj's, _fT and £ are functions of z only. Proposition 15.21 and Corollary 15.11 imply 
that these formulae hold only in the range C\N~^ < z < c 2 iV, where c\ and c 2 are constants 
independent of z and N. 

Equation (15. 15j) implies that 

C=^ + 0(M), A 2 ^0. (6.2) 

From this expression, the definition of £ (I4.2cp and the behaviour of the A^'s as v 2 — > (see 
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equation ( 15.22ajl ). we obtain 



8\0gGN Z 3,21 ^ / Z3 \ „ , „ r i ev 

Wz = 2 - ir ^ 3 + ( h ) + (2iV 5) ' (6 ' 3a) 

d log Gat 1 9 



log ((9 („(oo)+f- s ^3 7 )e(r*(oo)+ ? 



(9t 2 <9t \ V 2iriVN J \ y ' 2niVN 

+ l£-4)(i^(4)) + o(^S). ( a, b) 



Integrating these formulae gives 

r. / V 9 / 2 4 2e\ t 2 A^5 

= l^exp - - iV^3 _ ATT + 

\ 4 9 3 \ / 9.3 



5 4 



x (« („(oo) + c - ^) 9 («(oo) + c + 3^) ) ' (6-4) 



where B\ T is the constant 



B n = J rn (jl e'^^j Pgve(x u x N )d N x, (6.5) 

which is independent of t. 

The error term O yz^/Nij in equation ( 16. 4ft diverges unless z = O (iV 1 / 4 ). Therefore, 

we will restrict the validity of (16 ,4p to CiN~^ < z < C2N^. 

We now use the expression of the period II in equation (15.311) to simplify the theta 
functions in formula ( 16.41) . Combining the definition of the theta function (14.161) and the 
limit of u(oo) as A2 — ► ( 15.391) gives 

( U (oc) + c ± ^) = 1 + £ (^-f e 2 ™(j +t± ^ +0< ^») . (6.6) 
where cjq is a constant of order 0(1) in A2. Hence, we have 

* ( u(00) + < ± s&f) = 1 + (| eCJvA2 ) cos ( 2 - (i + < ± 5^37) ) + ( A3 J • ( 6 - 7 ) 

Substituting this formula back into (16. 4p and letting N — > 00 followed by e — > gives 

, , ^wOM) / 9 / 2 4 2e \ t 2 iVi\ /^ 2 \ . . 

lim lim ^ - ; exp ( N*z* - g-f = exp — . 6.8 

This limit completes the proof of Theorem 11.11 
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